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SUMMARY 


Bending-moment distributions were measured on a model helicopter 
rotor blade under hovering and simulated forward-flight conditions. 

The hinged-blade configuration was tested up to an advance ratio \i 
of 0 . 50 , whereas the fixed-at-root configuration was investigated up to 
and including \i = 0 . 90 - 

Curves of maximum-bending -moment distribution are presented for all 
test conditions. Harmonic bending moments have been found as a result 
of a harmonic analysis of the data for p = 0.22 and p = 0.^7. Theo- 
retical calculations have been carried out at advance ratios of 0.22 
and 0.47 Yor the hinged- and fixed-at-root conditions, respectively, 
and comparisons are made between experimental and theoretical results . 
Agreement between the results is reasonable in view of the assumptions 
made in the theory and the experimental errors involved. 

Resonance studies have been made on three sets of blades of differ- 
ent stiffnesses for the purpose of comparing the experimentally deter- 
mined resonance peaks with those indicated by theory. 

Aerodynamic -loading expressions are developed for the fixed-at-root 
blade which include the effect of elastic flapping. These results are 
used in the Goodyear method which is modified for the fixed-at-root 
condition. 


INTRODUCTION 

The purpose of the present program, which was conducted at M.I.T. 
under the sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics, was to investigate the possibility 
of measuring bending moments on a small-scale wind-tunnel model at advance 
ratios up to |i = 1.0. 
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Present theory permits the calculation of blade bending-moment 
distributions which are correct only to the degree of validity of the 
aerodynamic loading applied to the blade. At the present stage of 
development of rotor aerodynamic theory no method exists for the evalu- 
ation of the induced flow through the rotor in forward flight, and it 
is therefore necessary to approximate the aerodynamic loading. Conse- 
quently, the bending moment takes on a distribution and maximum value 
which is known to be in error. A reliable method of measuring bending- 
moment distributions on a rotor blade by means of wind-tunnel tests with 
small-scale models offers a possibility of checking present theory, 
determining the extent of the errors involved, and establishing the 
importance of higher harmonic bending moments not readily susceptible 
to theoretical analysis. , 

Since one of the important limitations on more general helicopter 
utility is the relatively low maximum speed attainable with flexible 
flapping blades, it may be desirable to extend the helicopter speed 
range possible by using rigid nonflapping blades. This immediately 
imposes a structural problem since the alternating stresses are expected 
to be high. This program has therefore also included tests beyond the 
present-day normal operating range for the purpose of contributing some 
information toward the evaluation of these stresses. 


SYMBOIS 

a slope of lift curve 

c blade-section chord, feet 

L rotor lift, pounds 

m mass of blade per foot of radius, slugs per foot 

r radial distance to blade element, feet 

R blade radius, feet 

Uip component at blade element of resultant velocity perpendicu- 

lar to blade-span axis and to axis of no feathering, feet 
per second ' 

Up component at blade element of resultant velocity perpendicu- 

lar both to blade-span axis and Urp, feet per second 

V induced inflow velocity at rotor, feet per second 


NACA TN 2626 


3 


a 

3 

G 


true airspeed of helicopter along flight path, feet per 
second 

rotor angle of attack, positive when shaft axis is pointing 
rearward 

blade flapping angle at particular azimuth position 
blade-section pitch angle 


tip-section ratio ( ^ 


mass density of air, slugs per cubic foot 
aerodynamic loading, pounds per foot 

blade azimuth angle measured from downwind position in direc- 
tion of rotation 


inflow ratio 


V sin g - V 
V 17R j 




rotor angular velocity, radians per second 


DESCRIPTION OF APPARATUS 


The model rotor under test in this program consisted of a three- 
blade system 5 feet in diameter. The blades had the NACA 0015 profile, 
a rectangular plan foim with no twist, and a chord of 3 inches. The 


rotor hub had horizontal hinges located at a radius of Ig inches to 

permit flapping motion and lag hinges with damping located at a radius 

of 2^ inches. Each blade had a 3/8- inch screw at its root coincident 

8 

with the quarter-chord line which fitted into a tapped split sleeve 
with clamping scfews on the hub which enabled the blade to be set and 
held at any desired pitch. 


The blade construction consisted of an aluminum tubular spar 
enclosed in a balsa profile as shown in figure 1. These pictures show 
the stages of assembly and the resulting aerodynamically clean surface 
of the instrumented blade. A 3/4- inch-diameter aluminum tube (l/32-in. 
wall thickness) was rolled down into one having an oval cross section 
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measuring O.3O by 1.00 inch. A root fitting was then inserted and 
riveted in one end and 10 pairs of strain gages were attached to the 
tube at radial stations shown in figure 2. Two balsa sheets were routed 
out so that one -half of the spar would fit into each piece. The three 
parts were then glued and clajnped together and the profile was cut by a 
special tool placed in a jointer. After sanding to the proper shape, 
the blade was balanced about the quarter chord by inserting small lead 
weights in the leading edge. The two dummy blades were constructed in 
the same manner but did not have strain gages attached to the spar. 

The resulting blades had a bending stiffness of l6,000 pound- inches 2. 

The slip-ring and brush assembly consisted of a total of seven 
rings and three sets of brushes. Six rings were made of coin silver, 
whereas the seventh consisted of Bakelite equipped with contacts for 
the purpose of indicating azimuth position. The silver graphite brushes 
were mounted in Bakelite holders and brush pressure was supplied by a 
coil spring fitted into a brass cap which screwed into the brush holder . 

In order to measure the bending moment at different stations suc- 
cessively without stopping the rotor between readings, a solenoid stepping 
switch was provided which made it possible to insert any desired pair 
of strain gages into the bridge. The axially symmetric switch was 
mounted coaxially with the center of rotation and the two leads from 
the solenoid were attached to two of the six silver rings. 

It was found that, in order to minimize the effects of slip-ring 
and brush resistance variations, it was necessary to mount all four 
arms of the resistance bridge on the rotor. Two d\immy strain gages were 
therefore mounted on small metallic blocks which in turn were mounted 
on inibber pads and placed close to the center of rotation. It should 
be noted that a pair of strain gages forming two arms of a Wheatstone 
bridge are mounted at each station so as to cancel out centrifugal 
strains and to measure pure bending. 

A Hathaway type MRC-12 strain-gage control unit was used for the 
purpose of amplifying the bridge unbalance. This unit is a compact, 
portable six-channel amplifier and power supply. The unit contains a 
2000-cycle-carrier oscillator element which supplies voltage to the 
bridge circuit. Its amplifier drives a Consolidated oscillograph galva- 
nometer type 7-112 having a sensitivity of 30 inches per milliampere 
from the voltage output of a strain-gage bridge. 

Since there were 10 strain-gage stations on the blade and only 
6 channels for amplification, it was necessary to have two stations 
alternately on one channel. A channel-selector control box was made 
with a stepping switch identical to the one on the rotor and synchronized 
with it so as to insert the gages into the proper channels . It should 
be noted that only one galvanometer was used since four slip rings 
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allow only one strain-gage "bridge to be energized at any instant. The 
switch in the channel-selector control box also connected the galva- 
nometer to the proper amplifier output channel. 

Two rotor mounts were provided for the hovering and tunnel tests. 
The conical hovering mount is shown in figure 3 where the rotor disk 
was approximately 10 feet above the floor. The drive motor was mounted 
with its shaft in a vertical position inside the cone and the rotor was 
connected to the motor by means of two universal couplings and a short 

shaft. Figures and 5 show the streamline mount for the 5- T^-foot 

wind-tunnel tests. Here the motor was slung underneath the mount and 
a comparatively long shaft was provided with a universal joint at each 
end. The mount could be tilted in the fore-and-aft direction to obtain 
the desired inclination of the shaft axis. 

PROCEDURE 


The sensitivity of the entire strain-gage and recording system was 
obtained by a calibration of the galvanometer deflection against known 
bending moments which were applied to the instrumented blade. The 
flapping hinge was locked so that the blade was horizontal and known 
forces were applied at the tip of the blade. The sensitivity of each 
station was thereby obtained in tenns of bending moment per inch of 
galvanometer deflection. The radial load was not simulated in general 
during calibration since centrifugal loads theoretically cancel at each 
station when the top and bottom strain gages occupy adjacent arms of a 
bridge. A check was made by putting a typical test blade in tension on 
a test rig with the result that there was not a complete cancellation 
of stresses probably because of small eccentricities. In any event if 
there were a possibility that the centrifugal loads would not cancel 
exactly, the effect would be taken into account in the zero reading 
recorded at each test rotational speed. 

The general procedure for obtaining bending-moment data was to set 
the blades at zero pitch and balance out 5 of the 10 stations successively. 
The rotor was then brought up to a desired speed and the galvanometer 
reading for each of the five stations was recorded. The same procedure 
was carried out for the remaining five stations. The pitch setting 
was then changed to that desired for the test and the procedure described 
above was repeated. In the tunnel tests, the zero pitch-setting readings 
were taken with the tunnel off so that only the effects of centrifugal 
force and blade weight were taken into account. The bending moment 
was then a measure of the distance from the line of the record established 
at 0=0^ and the trace at another pitch setting multiplied by the 
proper sensitivity factor. 
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With regard to blade twist due to the inclination of the blade 
principal axis of inertia to the plane of rotation, the effect was found 
to be small so no correction to 9 q was made. This was concluded from 
the knowledge of the blade torsional stiffness and a calculation of the 
torsional moments acting- on the blade at operating rotor speeds. 


DISCUSSION 
Inst rumen t at i on 


Since the hinged rotor blade described was expected to be subjected 
to a maximum bending moment of the order of only 10 inch-pounds, it was 
necessary to use a rather sensitive system in this investigation. A 
direct recording galvanometer (Consolidated type 7-112) in conjunction 
with an amplifier was found satisfactory and resulted in an over-all 
system sensitivity of approximately 10 inch-pounds per inch (ratio of 
bending moment applied to blade over oscillograph screen deflection). 

This figure varied slightly with the station and channel. Under these 
conditions the slip -ring and brush disturbances were at first intolerable 
but eventually were reduced, to a negligible level by insuring against 
dirt on the rings , providing a good running fit between brush and bushing, 
and eliminating any movement between the spring and cap in the brush 
holder by soldering the spring in the cap. When the fixed-at-root condi- 
tion was tested, a system sensitivity of approximately 40 inch-pounds 
per inch was used, since the moments were larger, with the result that 
the slip-ring hash was at an even lower level. 

A number of blades of different stiffnesses were tested in the 
hovering condition with strain gages at only one station inboard. A 
blade having a bending stiffness of 3OOO pound-inches^ was available 
whose construction consisted of an aluminum spar (1/8 by 3/^ in.) inserted 
in the leading edge of an NACA 0012 balsa profile. It was found that 
this blade undeiwent appreciable strain during rotational tests but very 
small bending moments resulted because of the low stiffness. The errors 
in the bending-moment measurements resulting from this condition were 
as much as 25 percent. Consequently, a blade having a stiffness of 

60.000 pound- inches^ was made in a manner similar to the final test 
blade described in the section "Description of Apparatus” using a steel 
tubular spar in an NACA 0015 profile with a pair of strain gages located 
at r/R = O.kO. The result was a blade with larger bending moments 
under operating conditions but smaller strains which were difficult to 
record with any degree of accuracy. From these considerations a blade 
having a stiffness of 20,000 pound-inches^ was decided upon and an 
aluminum-spar blade was constructed with a pair of strain gages located 
at r/R = 0.35* The actual stiffness on the resulting blade was 

16.000 pound-inches^ and it was found to be satisfactory in the hovering 
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tests since the degree of accuracy with which bending moments could be 
measured was within 10 percent. The final set of three blades was then 
made with this stiffness as discussed in the section "Description of 
Apparatus." Compared with the Sikorsky R-6 full-scale rotor blade which 
has a stiffness of 3«24 x 10^ pound- inches2 at the 50-percent- spar station, 
the model blade used in these tests is 9^ times stiffen based on the 
results of the development in appendix A. 


Test Results 

Hovering condition .- Typical hovering test results are shown in 
figures b(a) and 6(b) for the hinged- and f ixed-at-root conditions, 
respectively. Theoretical calculations have been carried out and the 
results are superimposed on the experimental results for the ease of 
comparison. A number of aerodynamic loadings were put into the Goodyear 
method (reference l) for calculating steady-state bending -moment dis- 
tributions of hinged rotor blades. Figure 6(a) shows the experimental 
curve and the theoretical distributions for uniform inflow with no tip 
loss, uniform inflow with 3-percent tip loss, and inflow variation with 
no tip loss. The effect of including a 3 -percent tip loss lowers the 
curve in general compared with the no-tip-loss curve and produces a region 
of negative bending moment in the neighborhood of the tip. Although no 
negative bending moments were obtained near the tip, the possibility of 
a bending -moment reversal occurring farther outboard than the last strain- 
gage station is allowed for as shown in figure 6(a). The inflow varia- 
tion (no tip loss included) was calculated according to Knight (refer- 
ence 2) and the modified aerodynamic loading was used in the bending- 
moment calculation. The resulting distribution conforms quite closely 
to that of the experimental result in the region of the outboard 50 per- 
cent of the blade but fails to agree with the maximum moment indicated 
by the tests. 

Figure 6(b) shows a comparison of theoretical and experimental rotor- 
blade bending -moment distributions for the hovering, fixed-at-root condi- 
tion. The theoretical calculations were made in accordance with the 
modified Goodyear method developed in appendix E of this report. The 
aerodynamic loadings used in this method were the same as for the case 
of the hinged blade. It can be seen from figure 6(b) that there is 
little difference between the results when using a tip loss of 3 percent 
and an inflow variation. The experimental distribution is shown to be 
in good agreement with the theoretical calculations. 

In the above calculations for the aerodynamic loading a two- 
dimensional lift-curve slope of 6.0 per radian was assumed for the hinged- 
and fixed-at-root conditions. 
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Simulated forvard flight .- The wind-tunnel tests for the hinged-^ 
at-root condition were run at a constant rotor speed, pitch setting, 
and shaft inclination. The tunnel speed was varied to obtain advance - 
ratios of 0.10, 0.22, O.3O, O.kO, and O.5O. It was found that a [x 
of 0.50 was the maximum that could be safely obtained because of an 
excessive amount of flapping and vibration. Figures 7(a) "to 7(a) are 
traces from the oscillograph records. The different scales on the 
ordinates of these tracings are due to the different sensitivities of 
the amplifier channels and bridges. For a given station the bending- 
moment trace varied within 10 percent from one cycle to the next for the 
hinged-at-root condition. An average cycle was therefore chosen for 
presentation in these plots. The pronounced fifth-harmonic bending 
moment is present because of. resonance with first-mode blade bending 
and contributed a large amount to the total bending moment. Stations 1 
to 3 ar^ not shown since the signals were small compared to the hash 
level. 

Figure 8 shows the family of maximum-bending -moment curves and the 
test conditions under which the data were obtained. These bending moments 
occurred at an azimuth angle of approximately 235^ for all the advance 
ratios tested. The bending moments decreased in going from a p of 0.10 
to p = 0.30 because of the rather fictitious conditions which were- 
set up in the tunnel since the effect of a variation in p alone (by 
changing tunnel speed and keeping rotational speed constant) on bending 
moment was desired. .The bending-moment traces show a decrease in the 
steady-state (mean) value as p is increased at a constant oc, 0, and 
rotor speed indicating a decrease in lift. However, as p is increased 
above O.3O, the first- and fifth-harmonic contributions have become 
larger and account for the increases in the magnitudes of the maximum- 
bending -moment distribution. 

The bending moment at a station on the blade may be expressed as 
a Fourier series: 

M = Mq + M^cos^ -I- Mti sin i ^^b sin 2>|^ + 

^3a 3^ ^3b sin 3t + . . . 

The coefficients of this harmonic series have been determined by a 
graphical method presented in appendix B. The results of a harmonic 
analysis for p = 0.22 are presented in figure 9* It can be seen here 
how large the fifth-harmonic sine coefficient is when compared with the 
other components. The second, third, and fourth harmonics were negligible 
and the fifth-harmonic cosine component is not shown since it was also 
small. It should be mentioned that more points do not appear on the 
plots in the vicinity of the tip since the moments are small here and a 
harmonic analysis is difficult and probably inaccurate under these 
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circumstances. It is felt that the points which determine the critical 
portion of the curve are presented. 

It was desirable to obtain some theoretical check on the data 
discussed above. Therefore calculations were carried out using the 
Goodyear tabular -dynamic method (reference l) to find the first -harmonic 
sine component of the bending -moment distribution. The steady-state 
component was also computed from this reference. Another method for 
calculating harmonic bending moments was investigated since it was con- 
sidered a relatively new approach to the problem which offered possible 
advantages over other methods. The De Guillenchmidt method (reference 3) 
is discussed in appendix C and calculation tables with an example are 
presented. The results of the above computations are presented in 
figures 10 and 11 where a comparison between the two theoretical approaches 
and the experimental results can be made. It should be emphasized that 
the same aerodynamic loadings were used in the Goodyear and De Guillenchmidt 
methods so that any differences in the results are due to the methods 
alone. In obtaining these aerodynamic loadings^ no tip loss,, constant 
induced flow, and a lift slope of 6.0 were assumed. The magnitudes of 
the corresponding bending -moment distributions are in reasonable agree- 
ment. However, the shapes of the theoretical and experimental cuives 
are different probably because of the aerodynamic loading which was 
assumed for the calculations. There is good agreement between theory 
and experiment in regard to the azimuth position at which the maximum 
bending moment occurs as noted in figure 10 . It may be seen from this 
figure that the magnitude of the experimental curve is larger than the 
theoretical result. The factor contributing to the discrepancy here is 
the fifth harmonic present in the experimental result but lacking in the 
theoretical result since no fifth-harmonic loadings were put into the 
calculations and therefore a fifth-harmonic bending moment could not 
possibly result. In general, comparison between theory and experiment 
indicates that present methods of computing blade bending moments are 
not satisfactory probably because of the inability of theory to provide 
the aerodynamic loading to make it possible to calculate the bending- 
moment distribution which actually exists on a given blade under operating 
conditions. 

It was also planned to carry out the wind-tunnel tests for the fixed- 
at-root condition at a constant rotor speed, pitch setting, and shaft- 
axis inclination, but as high advance ratios were approached it was 
found necessary to decrease the rotor speed to ^00 rpm in order to 
insure against overloading the blade. Data were therefore obtained at 
advance ratios of O.8O and O.9O at a lower rotor speed than the lower 
values of i-i, but, in spite of these precautions, on attempting to 
obtain \i = 1.0 the instrumented blade failed in fatigue and proceeded 
down the diffuser section of the tunnel. The aluminum-spar failure 
occurred just outboard of a steel insert at the root of the blade. 
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Figures 12 (a) to 12 (f) are the bending -moment traces from the 
oscillograph records for the fixed-at-root condition from p = O.15 
to u = 0.90. The scale of the vertical axis (bending moment) varies 
from one station to another again because of the difference in the sensi- 
tivities of the amplifier channels. The cycles presented in these fig- 
ures are considered typical of those occurring over the short time inter- 
val during which records were made. It should be noted that the cycles 
for a given condition repeated with negligible difference for the fixed- 
at-root condition. Data for stations near the tip are not shown since 
the oscillograph records showed mere ripples, indicating very small 
bending moments. At advance ratios of O.8O and O.9O the signal from 
station 10 was not available since the strain gage at this station failed. 

Plots from the oscillograph records are given in figure 13 which 
shows the variation of bending moment at station 9 with azimuth position 
for the values of advance ratios tested. It may be noted that the 
azimuth position at which the maximum moment occurs changes as p 
increases. The drop in bending moment from p = O.60 to p = O.80 is 
in agreement with the decrease in rotor speed and hence aerodynamic loading. 
However, at p = O.9O with the lower rotor speed, the maximum bending 
moment exceeds that at p = O.60 because of the increase in tunnel 
speed which resulted in an increase in aerodynamic loading. It is sus- 
pected that the maximum bending moment on the blade at failure was in 
the neighborhood of 200 inch-pounds which corresponds to a tensile 
stress of 20,000 poiinds per square inch. 

The oscillograph records for all the runs were analyzed for the 
pxirpose of obtaining the maximum spanwise bending -moment distributions 
and the azimuth positions at which they occurred. The results of this 
analysis are shown in figure lU where the maximum spanwise bending-moment 
distributions for various p conditions are given. Implicitly, these 
curves show the change in spanwise moment distribution for corresponding 
changes in tunnel or simulated forward- flight velocity for the fixed-at- 
root rotor configuration up to a p of O.60. The curves have been 
extrapolated to the center of rotation since no strain gages were located 
on the hub. The curve for p = O.8O is below that for p = 0.60 
because of the particular combinations of tunnel and rotor speeds used 
to obtain these advance ratios. 

A harmonic analysis was carried out on the bending -moment traces 
for the p = O.i^-7 condition and the results are shown in figure 15. 

It is of interest to note that the first-harmonic components are of the 
same order of magnitude but of opposite sign and that there is a large 
contribution from the second-harmonic cosine component. The higher 
harmonics were foiind to be negligible. 

In order to check the experimental results against theory for the 
fixed-at-root condition, it was first necessary to derive some convenient 
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expressions for the aerodynamic loadings under these conditions. This 
development is given in appendix D and considers only first -harmonic 
effects. The Goodyear method (reference l) has been modified in 
appendix E for the purpose of calculating steady-state and first-hai^nonic 
bending -moment distributions for the fixed-at-root blade. Also included 
in appendix E is a method of iteration which is used to cope with the 
problem of coupling between the aerodynamic -loading components. 

A comparison of the theoretical and experimental harmonic -bending- 
moment distributions for the fixed-at-root condition at p = 0.47 is 
presented in figure 16. There are a number of factors here which lead 
to the discrepancy between theory and experiment, namely the assumed 
constant induced velocity, the neglect of the reversed- flow region, and 
a probably high two-dimensional lift-curve slope (a = 6.0). In addition, 
no wind-tunnel wall corrections were applied since there are no reliable 
methods available by which such corrections can be estimated, particularly 
as they may affect harmonic bending moments . A comparison between 
maximum bending moments as discussed for the hinged-at-root case was 
not made since the second-harmonic cosine component of bending moment, 
which apparently plays a large part in determining the total bending 
moment on the fixed-at-root rotor blade, could not be readily calculated 
with the present development. 


Resonance Studies 

The problem of rotor-blade flapping vibration in the hovering 
condition was met with early in the investigation and an attempt was 
made to understand the problem and alleviate this condition. Rotor 
balance and shaft alinement were Improved with the result that blade 
vibration was reduced only a small amount. It was noted that the 
vibration was a function of pitch setting, the effect being rather small 
at zero pitch. This indicated the source of trouble to be of an aerodynamic 
nature and further investigation showed that the vibratory frequency 
(cycles per revolution) of the rotor blade varied with rotational speed. 

A resonance condition was suspected. 

A theoretical sui^ey of blades having three different stiffnesses 
was made by carrying out a calculation of their natural bending frequencies 
as a function of rotational speed by the method of reference 4. The 
results of this calculation for the hinged- and fixed-at-root conditions 
are shown in figures l?(a) and 17(b), respectively, for blades having 
stiffnesses of 3OOO, l6,000, and 60,000 pound-inches^. The intersections 
of these curves and the straight "per rev" lines are resonance conditions. 

The resonance problem was further investigated to obtain a better 
understanding of the phenomenon. Tests were run on three sets of blades 
having stiffnesses of 3000, l6,000 and 60,000 pound-inches^ in the 
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hovering condition. Oscillograph records, taken at small increments 
in rotational speed, gave the strain-gage response to blade vibration 
at pitch settings of 0 ^ and 12 ^. These records vere analyzed to obtain 
the steady-state or mean bending moment and the vibratory or resonance 
bending moment at a given rotor speed. Figure 18 gives two sample oscil- 
lograph records showing the resonance phenomenon and the above terms 
defined. It should be noted that the vibratory motion of the blade was 
not of constant magnitude but consisted of an irregular increase and 
decrease of the double amplitude at* a constant rotor speed. In each 
case the maximum double amplitude has been considered and plots of the 
ratio of the maximum double amplitude of resonance bending moment to 
steady- state bending moment against rotational speed for different blades 
are shown in figure 19 for the hinged-at-root condition and figure 20 
for the fixed-at-root condition (p = 0 °). Superimposed on each plot 
is the appropriate natural- frequency spectrum for convenience in com- 
paring theory and experiment. Each peak of the resonance curve is 
labeled with the cycles per revolution obtained as a result of the oscil- 
lograph record analysis. This quantity refers not only to the peak but 
also to the neighborhood of the peak. For example figure 19 (b) indicates 
an experimentally detennined peak resonance condition of 6 cycles per 
revolution at 58O rpm, of 5 cycles per revolution at 735 rpm, of 4 cycles 
per revolution at 1000 rpm, and of 3 cycles per revolution at 1350 rpm. 

The intersections of the ”per rev” lines and the theoretically determined 
first -natural-bending-frequency curve occur at rotational speeds very 
close to the above except in the case of the resonance condition of 
3 cycles per revolution. This agreement carries through for the hinged- 
at -root-condition tests (fig. 19)* 

The calculations of the first natural bending frequencies for the 
fixed-at-root 'rotor -blade condition were based on a blade length of 
2 feet instead of 2.5 feet since the root fitting had a stiffness of 
approximately 600,000 pound-inches^ and the actual blade root was 
located at a radius of 6 inches . The agreement between theory and 
experiment is good for the low-stiffness blade but becomes worse as 
blade stiffness is increased. This root condition introduced contri- 
butions to the effective stiffness of the rotor system from the root 
fitting and flexibility of supporting structure which were not considered 
in the theory. It seems reasonable to assume that these effects would 
be of greater importance as blade stiffness is increased and this appears 
to be borne out by experimental results presented in figures 20(b) 
and 20(c) . 

A point of interest to be noted from figures 17 to 20 is the vari- 
ation of the ratio of resonance bending moment to steady- state bending 
moment with blade stiffness. For the hinged-at-root condition, the blades 
having stiffnesses of 3000, l6,000, and 60,000 pound- inches^ indicate 
maximum resonance-peak ratios of 1 . 0 , O.38, and O.60, respectively. 
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Also for the same blades fixed at the root, the maximum resonance-peak 
ratios are 0.37, 0.20, and 0.^7, respectively. In each condition note 
that the blade of relative medium stiffness experiences the lowest 
resonance level. It appears from this brief study that there is an 
optimum stiffness which will minimize blade resonance; however, it is 
felt that a general statement of this nature cannot be made at this time 
without a more thorough investigation of the resonance phenomenon in 
helicopter rotor blades . 


CONCLUSIONS 


Bending-moment distributions were measured on a model rotor blade 
under hovering and simulated forward- flight conditions in both the 
hinged- and fixed-at-root configurations and a comparison was made with 
theoretical calculations. It was found that: 

1. For the hovering, hinged-at-root condition, the best agreement 
between theoretical and experimental bending-moment distributions was 
obtained when the inflow ratio variation according to NACA TN 626 was 
considered. No measurable amount of negative bending moment at the tip 
was found which indicates less tip loss than is usually assumed. 

2. In the hovering, fixed-at-root condition, the experimental bending- 
moment distribution fell within the region determined by theoretical 
calculations obtained as a result of inserting different types of aero- 
dynamic loading into the bending-moment calculations. 

3. The fifth-harmonic bending moment due to resonance in the forward- 
flight, hinged-at-root condition is relatively large and accounted for 
much of the discrepancy between the theoretical and experimental total 
maximum bending moments . An appreciable difference between the general 
shape of the distributions predicted by theory and measured experimentally 
exists and is probably due to the distribution of aerodynamic loading 
assiamed in the currently available theory. 

k. In the case of the fixed-at-root, forward-flight condition, it 
was found that the second -harmonic bending moment was the largest com- 
ponent contributing to the total moment. The agreement between theory 
and experiment in regard to the steady-state and first-harmonic bending 
moments for an advance ratio of 0.47 vas reasonable in view of experi- 
mental errors and the assumptions made in the theory. 

5. A brief study of rotor-blade resonance phenomenon was made in 
the hovering condition. Good agreement between theory -and experiment was 
obtained for the hinged-at-root condition since it was found that the 
experimental resonance peaks occurred at rotational speeds very close 
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to those where the natural-hending- frequency curves crossed the lines 
of cycles per revolution. Agreement in the case of the f ixed-at-root 
condition was reasonable in view of the contributions to effective 
stiffness of the rotor system from the supporting structure. 

6 . Information regarding the actual aerodynamic loading or induced 
flow and blade resonance phenomenon is needed if the results of bending- 
moment calculations are to represent the actual situation existing on a 
rotor blade under operating conditions. 


Massachusetts Institute of Technology 
Cambridge, Mass., April 3^ 1951 
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APPENDIX A 

BASIS FOR COMPARISON OF STIFFNESSES OF VARIOUS ROTOR BLADES 


Consider a blade under the loading W = f(r). The deflection S at 
the tip of a cantilever beam is; 


6 


_ 

■ 3 EI 


where R is the blade radius. Now 


W(r) = ■^peuXj.cfi^r^ dr - m' dr fi^r 
Let X = — so r = xR and dr = R dx. Then 


^ ^ ao^cfi^R^x^ dx - m*^^^R^x dx 

2 

Now is mass per unit length and oc L where L is the 

dimension of length. Also c and R have dimensions of length L. 
So 

W a 


From the deflection expression there is obtained: 

& ^ cc 

R El El 

Equating the expressions for 6/R of the model and the full-scale 
blade, there is obtained: 


(2R)/l^^ _ (nR)p%^ 

(ei)m (ei)f 


l6 
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where the subscripts M and F denote model and full scale, 
respectively. Hence, 




(El)j, 


(r2R)p%^ 


Comparing the R-6 "blade and the model hlade under test there is 
obtained: 


h? = 

Lp 19.0 


= 0.1315 


— = 2.98 X 10"^ 


.1^ . 0.176 
i^^)J ( 500 )^^ 

F 


From reference 1 , 


P 

(EI)„ = 3-24 X 10 lb-in. 

Jb 


at X = O.5O; hence, 

(EI) = 3.24 X 10 ^ X 0.176 X 2.98 X 10"^ = 170 lb-in. ^ 

M 

This result means that a model of the R -6 rotor blade having a diameter 
of 5 feet and rotating at 8OO rpm should have a stiffness EI of 170 poimd- 

inches^. The EI of the model blade Investigated in this report was 

16,000 pound-inches^ and is therefore 94 times as stiff as a model of 
the R-6 when rotating at the above rotational speed. 
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APPENDIX B 

GRAPHICAL METHOD FOR HARMONIC ANALYSIS 
OF A PERIODIC FUNCTION 

Method.- Let the c\arve such as the one shown in the accompanying 
figure be represented by the following series: 

y = bQ + a^^ sin x + b^ cos x + E 2 sin 2x + b 2 cos 2x + . . . + 

aj, sin nx + bjj cos nx + . . . 



The curve is divided into 12 equal parts such that the distance 
between ordinates h is «/6, and the same ordinates are used in 
finding all harmonics. 

In order to find the first-harmonic coefficients a-^ and b-^, the 
following construction is made: Starting from a point 0, — y^ is 

laid out in the appropriate direction on the horizontal; y^^ is added 
to i y^ at an angle of 30° to the horizontal; yg is added at an 

angle of 60°; and so forth. Note that when the ordinate is negative, 
it must be laid off in the opposite direction as shown in the fig\ire 
on the following page. 
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2 ^12 Jl 



The coefficients and h^ are found from the relations: 


®1 ■ g 


= I 

To find the coefficients of the second harmonic a similar construe- 
tion is made, but the ordinates are laid out at intervals of 60° instead 


of 30 • The expressions for ag and b2 


^2 ■ 6 ^ 2^2 


b2 = g ON 2 


are: 
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In general, the nth harmonic may he found by construction where the 
ordinates of the curve are laid out at intervals of n(2jt/M) where M 
is the number of parts into which the period of the curve is divided. 
The expressions for the coefficients are: 


^ M ^n^n 


^n ~ M 

Proof of above construction .- The trapezoidal rule states that the 
total area \inder a curve is: 




+ y + y + 
J-1 J'g 


where h is the distance between ordinates, or h = 2 jt/m. 
Now the Fourier series coefficients can be written: 


^1 - n. 


y sin X dx 


For the integral substitute the expression for the area given above: 


= -^ yQ sin 0° + y^ sin 30° + y^ sin 6o° + 


■ 2 ^12 


or 


1 ^ Ifi ^0 ^1 ^2 ’ 


+ i ^12 360°] 


Also 


b 


1 



cos X dx 
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So 

"" i(i ^0 30° + 72 cos 60° + . . . + -| 73_2 cos 360° j 


Now the graphical process described above is that indicated in the 
parentheses, thus: 


^1 - M 


bi = ^ ONi 


Also 

^0 ^ ^ "'in^ ir(i ^0 ^1 ^2 i ^ 12 ) 

Hence 


_ 1/1 

^0 M\2 


yo + 7i + 72 



NACA TN 2626 


21 


APPENDIX C 


DISCUSSION OF DE GUILLENCHMIDT METHOD FOR CALCULATION 


OF BENDING STRESSES IN HELICOPTER ROTOR BLADES 


The development of this method for calculating bending stresses in 
rotor blades is given in reference 3. The procedure for the general 
method, the method as applied to the model rotor tested in this investi- 
gation, and sample calculations are discussed in this appendix. 


The procedure for calculating the elastic bending moment is as 
follows : 

1 . Calculate the natural frequencies and corresponding mode shapes 
of the nonrotating blade. Reference 3 states that the rotational and 
nonrotational natural modes are the same within a few percent. A cal- 
culation has been carried out on the R— 6 rotor blade and the results are 
presented in figure 21 where it can be seen that the first natural mode 
does not vary greatly with rotational speed. The frequencies and mode 
for the model rotor blade were obtained by the use of the Myklestad 
method (reference 5, PP- 184-203). The rotating natural bending fre- 
quencies were calculated from the results of reference 4. 

2 . Calculate the aerodynamic loading (reference l), the centrifugal 
loading, and the total external loading on the rigid blade at different 
azimuth positions. For the example illustrated here, eight positions 
spaced 45*^ apart were used as shown in tables I, II, and III. Curves of 
total external loading F^' are shown in figure 22 . No tip loss was 
included in the aerodynamic loading. 


3. Then evaluate the integral / dx at the corresponding 

'^0 

azimuth positions as illustrated in table IV. The mode functions t]]_ 
and t \2 for this example are shown in figure 23. 


Procedure for Calculation of Elastic Bending Moment 


4. The integral 



0 


then evaluated as shown in table V. 
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5. Thus it is possible to calculate (table VI) for the eight 
azimuth positions. The plots of g]^ and g2 are shown in figure 24. 

A harmonic analysis is performed on the g function using the graphical 
method described in appendix B of this report. This determines the 


of the Myklestad computation since this integral is the ratio of the 
rigid -blade bending moment to the square of the natural frequency. This 
distribution and the elastic -blade bending -moment distribution broken 
down into harmonics are shown in table VII. 

7 . The elastic bending moment obtained from a consideration of the 
second natural mode of the blade is shown in table VIII. The bending 
moment resulting from both modes is shown in table IX. It was found 


The De Guillenchmidt method is particularly useful for the investi- 
gation of the effect of a series of loading conditions on a given blade . 
This is true because^ once the natural frequencies and modes are deter- 
mined^ relatively little calculation is necessary to obtain both the 
harmonic- and total -bending -moment distributions for a given loading 
condition. 

For the blade tested in this investigation, the calculations showed 
that the effect of including the second natural mode was negligible in 
the case of the steady-state component, whereas the effect on the sine 
and cosine components was appreciable. However, for full scale the 
contribution of the various modes to the bending moment must be investi- 
gated for each blade design. 


values 




is one of the byproducts 



Comments on De Guillenchmidt Method 
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APPENDIX D 

CALCULATION OF EXIEESSIONS FOR AERODYNAMIC LOADING 
ON A FIXED-AT-ROOT ROTOR BLADE 


The rotor blade discussed in this appendix is one which is fixed 
at the root but is allowed to bend under applied loads because of its 
elastic properties. 

The lift on a blade element is: 

AL = ipac + UpUp^Ar (D1) 

where 


Urp = 4 - mE sin Q.t) (D2) 

and 

Up = X,® - z - pOR ^ cos m; (D3) 

^ dr 

where z is the deflection of a blade element from the plane perpendicular 
to the axis of rotation and passing throijgh the blade root. Now z can 
be written as: 


z = r(3 + y-j + sin fit + Yjjj cos Gt + . . . (d4) 

Only the first harmonic will be considered in this development; 
therefore: 


2 = yijfi cos at - yiiifi sin nt 


(D5) 


and 


dz 

dr 


^Yl *^11 *^111 


COS Qt 


dr 


dr 


dr 


(D6) 


2h 
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Hence 


Up = cos fit + yjjjfi sin fit - 


I 

HfiR cos fitm + + 


^11 *^111 \ 

— — sin fit + — T— cos fit 
dr dr / 


(D7) 


Substituting these quantities into the lift expression, equation (D1) 


T = 


Ar 


^ac-j 


0fi^(r + sin fit)^ + 


“ ^11^ cos rrt; ^xii^ sin ~ 


/ <3yj *^TT III 

nfiR cos fit 3 + -rr^ + sin fit + cos fit) 




dr dr 


dr 


(fir + M-fiR sin fit) 


Upon expanding and factoring out of there is obtained: 
T = ^acfi^^er^ + 20^rR sin fit + e^i%^sin^fit + XrR + 


Xp.R'^sin fit - yjjr cos fit - sin fit cos fit + 

y^^^r sin fit + yJJJ^® sin^fit - |irR3 cos fit - n%^3 cos fit sin fit 


dyi p p dyj “^yil 

lirR — — cos fit - M- R — — sin fit cos fit - M-rR — r — sin fit cos fit - 
dr dr dr 


p P *^^11 P p 

|i^R^ cos fit sin fit - M-rR — 77;:— cos^fit 


dr 


dr 




dr 


sin fit cos fit 


(D8) 
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This can. be written as follows: 


T = ^acf2^ \6r^ + 20|arR sin Q.t + % cos + 


Ytt 

X,rR + XkiR^ sin ftt - y^-jr cos fit ^R sin 2fit + Yuj^ sin fit 


(IR - cos 2Qt - ^3rR cos fit - ^ i-i^pR^sin 2fit - 


p.rR 


cos at - R^ sin 20t - li rR sin 20t - 

dr 2 dr 2 dr 


u *^^111 

cos at + ^ cos 30t - 2 I’R “dT“ 


r 2 x>2 

T~ dr 


<iy 


III 


I ~dF~ 20t - n- 


. r2 5^111 3in at - r2 ^ 


sin 3^^) 


Comparing equation (D9) with: 


T = Tj + sin at + cos at + 


the following expressions are obtained: 
The steady-state loading: 


Ti = •|paca2 


0 r^ + |j2r2 + 


XrR + 2 (^111 


<iy 


- r 


III 


dr 


The first -harmonic sine component: 


1 2 1 ^ o2 ^III] 

Tjj = •^aca l20rpR + XpR + “ Ij~ ,32. j 


(D9) 


(DIO) 


(Dll) 


(D12) 
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The first -harmonic cosine component: 


Tj-j-j — "^pahn 


-yjjr - 


dyi£ 


dr 


(D13) 


The above expressions are used in the application of the method of 
calculating rigid-blade bending moments developed in appendix E. 
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APPENDIX E 

MODIFICATION OF GOODYEAR METHOD FOR CALCULATING 
FIXED-AT-ROOT ROTOR-BLADE BENDING MOMENTS 
DiBCUBSlon of Modification 


In tills appendix the Goodyear method (reference l) for calculating 
hinged-rotor-blade bending moments will be modified so that rigid-blade 
moments may be determined. 

First consider the steady-state analysis for a rigid blade (i.e., 
fixed at root but allowed to bend) . The total shear at any point on the 
blade is: 



(El) 


where 

6 aerodynamic shear at a given point 

F centrifugal shear at same point 

M bending moment at given point 

z deflection from plane perpendicular to axis of rotation and 

passing through blade root 



Now 


z = rP + y^ + y^j sin Qt + Yjjj cos Dt + . . 
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So 


dz 

dr 



^yii 

dr 


sin fit + 


dr 


cos (2t . 


where y is the deflection measured from the 

The steady-state component of dz/dr is 
tion (El) becomes: 


undeflected position. 

3 + Hence eqiia- 

dr 




+ S 


(E 2 ) 


This equation can be approximated by dividing the blade into an 
Integral niimber of parts each of length X. Then dM/dr at the midpoint 
of an interval is: 


(m\ ^ ^i+1 ~ ^i 

W/i+l" X 


Hence the above equation becomes: 


M. 


i +1 


= M, + F X (p + + XS , 

^ i4V 


Let 


Mj^ = aj_ + bj^Z 


where Z = Mq = Bending moment at blade root. 
At the blade tip: 


^ = ^n + ^n^ = 0 


or 



(E3) 


(e 4 ) 


(E5) 
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Hence the bending moment at the root will be known once the 
values a^^ and b^^ are determined. To find the expressions for 

the a's and b's proceed as follows: 

At the blade root: 


and 


Mq - aQ + b^Z 


Mo = Z 


Comparing these two equations there is obtained: 


ao = 0 


bo = 1 


Now 




Then equation (E 3 ) can be written: 


%+i = ^i^ i 


Also equation (eU) can be written: 


/ a.X^ ^ b.X^z\ 

F + 2- vf- Z_ 1 


%+l ~ ^i+1 ^i+1^ 

Equating coefficients of equations (E 9 ) and (ElO): 


^i+1 = aj^ + XS 1 ^ -|\P^ 


y 

Till 


i+| i+|V J=o 


(e6) 


(E7) 


(e8) 


(E9) 


(ElO) 


(Ell) 
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\ 

'i+1 - ^ 1 Z_ EI7 

2-3=0 .3 


(EI2) 


Now consider the first -harmonic bending-moment analysis. The shear 
equation now reads: 


.s + Jp + n2 r Dwdr 


(E 13 ) 


where Jq is the inertia shear at the root and the quantity 

2 r 

Jq + I my dr is the inertia shear at r. The subscripts denoting 

Jo 

the harmonic have been dropped since the solution takes the same form for 
the sine as for the cosine component. 


Now j my dr can be written as: 

/o 

i 


I 

j=i 




(ei 4) 


and 


dr 


dz\ ^ ^ 


II 


dr 


(E15) 


So 


M 


i+1 




1+2 ^ ^+i i+2 J=I 


± 

I 


2 2 

X n m^yj 


(ei6 ) 


Let 


Mj^ = at + L^Mq + c^0 


(EI 7 ) 
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and 


+ ViMo + W^0 


where 0 = XJq. Now 


2 

X M, 




= 


El. 


0 . 

^EIq EI^ 




i-1 X^M. 


El. 


j=l 


So 


X% ^ 

and equation (EI8) gives: 

+ V^Mq + w^0 


By comparison 


TT X 

^1 = eT 


^l=Ei ^0 


„ X^ 

% = eT ""0 


(ei8) 


(E19) 


Proceeding further, equation (E19) gives: 


or 


y2 = 


El 


El 


xV 


El 


d d 

^2 = ii(ao + ^0 % + Cq 0) + |j[aQ + a^ + (b^ + b^M^ + (cq + cj^ 
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Wow 


By comparison 


72 = U 2 + V 2 M 0 + W20 


Uo = 


- % i?(^o ^3) 


''2 = I'l + ir(t’o + »x) 


«£ = «i * ir(=o * =i) 


It can be concluded that in general: 


U 


1+1 




^i+l = Vi 


^1+1 = + 


J— 1 


X^b. 




X^c. 


J=1 


(E20) 


Expressions for the a's, b's, and c's are found in the following 
manner: 

Consider the bending moment at the blade root: 


Mo = Mo 


Equation (EI7) becomes 


Mo = aQ + bQMQ + Cq 0 
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So by comparison 


ao = 0 


^0 = ^ 
Co = 0 


Now consider equation (E16) with i = 0: 


“1 - "0 * " ^ 1/2 * 


Since 


then 


W/1/2 


x^iyu 

«1 = «0 + Fi/ 2 + XSi/2 + t 


Comparing this equation vlth 


M2^ = 4- h]_MQ 4- Cj0 


ai - XSj^yg 

»1 - 1 * fi/2 K 


there is obtained: 
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Repeating this process for i = 1: 


2 

Mg = Ml + |j.^Mq + Ml) + + 0 + + ViMq + W 


or 


2 

Mg = ai + biMQ + Ci 0 + F 1 ^{mq + a^_ + biMQ + Ci^) + XS 1 

I2 Ig 

x^nS(^i ^i“o ^ ^1*^) 


Since 


then 


Mg = Eg + bgMQ + Cg0 


^2 “ ^*1 ^ 1 ^ 1 


X^ai 


El 

2 2 


X^^lUi 


bg = bi + F^l Ij.^ 1 + bi) + 


^ 2 ^ ^l^\l |[(Cl) + 1 + X^AiWi 


It can be concluded that in general: 


H+1 


H+1 


‘ "1 * ‘i+i ^ ii " 

2 j=i j 2 j=i 

= b, + F 1 ^ + X^tl^ ^ m.V, 


H*i ‘ '=1 * f lY 1 * XV X -j«, 

2 «J J=1 


I**) 


0 + 


(E21) 
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It is knovn that 

J=n 

XJq = - = 0 


as defined earlier where n is the tip station. Substituting equa- 
tion (E18) for in the above: 


XI 

0 = -II 


j=i 


+ M, 




__ xSn^m.V, +0 

j=l ^ ^ J=1 




(E22) 


The summations are determined from the calculation. Let them he: 


= Lx2fi2m^Uj 
N2 = IX^n^jVj 

"3 = 


(E23) 


Then equation (E22) becomes: 

0 = -Ni - N^Mq - N30 


or 


0(1 + N 3 ) + W^Mq = -Ni 


(E24) 


Also^ at the tip. 


^n " ^n ^n^ ® 


Hence 


0 = 


- - — Mq 
n n "" 


Substituting this expression for 0 into equation (E2U): 

(1 + '< 3 )(- ^ « o ) + " 2«0 = -»1 


(E25) 
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Therefore 



(E26) 


Knowing Mq^ 0 may be obtained and hence from the expression 

^i “ + ^i^ 


Discussion of Sample Calculations 

Steady-state component . - Table X, which shows the calculation for 
this condition, is very similar to that given in reference 1. Column (5) 
is calcxxlated in the same manner as in reference 1. The equation for 
finding the aerodynamic loading for the rigid blade is given in appendix D 
of this report. Likewise, column (7) is computed as in reference 1 . 

Notice that Ym must be neglected in the calculation of tj 
since y^^^ is not known at this stage of the computation. It was 

found that after completing a cycle (i.e., sine and cosine components), 
the effect of Xjjj °u was small. Therefore the result obtained 

from table X was taken as the steady-state bending-moment distribution. 

The computation check (table XI) for this calculation is identical 
to that given in reference 1 , and is the same for the first -harmonic 
components. 

First -harmonic components . - Table XII shows the calculation of a, 
b, and c. It should be noted that only a changes with the aero- 
dynamic shear, b and c remaining the same for given blade constants 
and rotational speed. 

The first-harmonic sine component was first calciilated by using the 
expression for Tjj given in appendix D neglecting Xjjj since it was 

not known at this stage of the calculation. The computation checi 
yields in addition the deflection yjj. 

Having yj and Yj^y it is then possible to calciilate Tj-[-j and 

the corresponding aerodynamic shear. A calculation of the a's according 
to table Xll(a) permits the computation of the cosine component of bending 
moment and consequently Yjjj in the computation check. 
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The sine component of the aerodynamic loading is then 

recalculated using y^-^j- just obtained. For this example it was found 
that yiii had a considerable effect upon Tjj and consequently on the 
bending moment. The iteration process is continued to find a new y^j: 

To decrease the manber of cycles in the iteration process the mean value 
of Yjj resulting from the first and second cycles is found. This 

value of yjj is then used in the second-cycle calculation of ’’’jjj* 

The res\iltlng cosine bending moment and consequently the distribution of 
yill are then found. This makes it possible to recalculate the sine 
component and find y^j which is compared with the mean value assumed 
above. In this example the two values of yjj were within 10 percent 
so that iteration was not continued. 

A plot of the bending-moment distributions during the above process 
for the blade under test is shown in figure 25. It shoixld be stressed 
that only the a's need be recalculated in carrying out the iteration 
since b and c do not change with aerodynamic loading. 

It has recently come to the attention of the author that the Cornell 
Aeronautical Laboratory has extended the Goodyear tabular-dynamic method 
to take into accoimt the additional aerodynamic loading introduced by the 
elastic deformation of a hinged rotor blade. This additional effect is 
considered immediately in the moment equation and the sine and cosine 
components are carried along simultaneously so that no iteration is 
necessary. The application of this method^ with the appropriate modifica- 
tions, to the fixed at -root blade has not been studied in detail, but it 
appears that such an application is possible. For the present case, 
where no possibility of resonance with the assumed loadings existed, the 
iteration -process appears to be a reasonably rapid and accurate method of 
estimating the bending -moment distribution on a semirigid nonflapping 
blade . 
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TABLE I 


VALUES OF AERODYNAMIC LOADING F^' (LB/FT) 
AT EIGHT AZIMUTH POSITIONS 


[calculation of F^' may be set up the same way as that given in reference iT] 


Station 

0 

0 

II 

i = 45° 

0 

II 

i = 135° 

ijr = 180 ° 

0 

IT\ 

OJ 

OJ 

II 

0 

0 

OJ 

II 

y = 315° 

0 

1 

0.098 

-0.136 

0.182 

-0.014 

0.270 

0.504 

0.550 

0.382 

1 

.138 

.313 

.496 

.581 

.516 

.341 

.160 

.073 

2 









3 

00 

VJl 

0 

1.27 

1.57 

1.58 

1.29 

.870 

.571 

.560 

2.25 

2.66 

2.93 

2.91 

2.61 

2.20 

1.93 

1.95 

k 

U.34 

4.52 

4.63 

4.60 

4.46 

4.28 

4.17 

4.20 

5 

7.11 

6.85 

6.66 

6.66 

6.85 

7.11 

7.29 

7.30 

6 
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TABLE II 

CALCULATION OF CENTRIFUGAL LOADING Fc' (LB/fT) 
AT EIGHT AZIMUTH POSITIONS 
[Fc' = Fc3] 


Station 

Fc 

(Ib/ft) 

i = 0° 

p = 0.0007 

= 45° 

p = 0.0056 

y = 90° 

p = 0.0444 

y = 135° 
p = 0.0930 

i|f = 180° 
p = 0.123 

t = 225° 
p = 0.117 

if = 270° 
p = 0.0778 

t = 315° 
p = 0.0293 

0 

8.87 

-0.0062 

0.0497 

0.394 

0.825 

1.09 

1.03 

0.690 

0.259 

i 

26.6 

-.0186 

.149 

1.18 

2.47 

3.27 

3.10 

2.07 

.777 

2 

0 

kk,2 

-.0309 

.248 

1.96 

4.11 

5.43 

5.15 

3.44 

1.29 

j 

62.0 

-.0A3A 

.347 

2.75 

5.77 

7.62 

7.23 

4.82 

1.81 

k 











79 . T 

-.0558 

.446 

3.54 

7.41 

9.80 

9.29 

6.20 

2.33 

5 

97.5 

-.0683 

.546 

4.33 

9.07 

11.9 

11.4 

7.59 

2.85 

6 
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TABLE III 

CALCULATION OF TOTAL EXTERNAL LOADING Fd' (LB/FT) AT 
EIGHT AZIMUTH POSITIONS 


[Fd’ = Fa’ - Fc] 



I 


h2 


NACA TN 2626 


TABLE IV 


CALCULATION OF 



Fd’Tii <ix 


AT EIGHT AZIMUTH POSITIONS 


Station 

’ll 

^|r 

= 0° 

^ = 


Mr ^ 

= 90° 

Mr = 

= 135° 

Fd'ni 

Fd’Tii X 5 " 

Fd'm 

Fd’Hi X 5 " 

Fd'Hl 

Fd'Ri X 5 " 

Fd'Hi 

^d '^1 5 " 



0 

-0.210 

- 0.022 

- 0.110 

0.0391 

0.196 

0.0445 

0.223 

0.170 

0.85 

1 

-.560 

-.088 

-.440 

-.0913 

-.'»57 

.384 

1.920 

i 

1.06 

5.30 

2 










3 

-.680 

-.599 

-2.99 

-.694 

- 3.‘^7 

.266 

1.33 

1.72 

8.60 

-.U62 

-1.058 

-5.29 

-1.067 

- 5 . 3 »^ 

-.082 

-.41 

1.32 

6.60 

k 










c 

.010 

.Ol^l^ 

.220 

. 04 l 

.205 

.011 

.055 

-.028 

-.ito 


.610 

It. 379 

21.91 

3.84 

19.20 

1.42 

7.10 

-1.47 

-7.35 

6 


Jo 

lx = Z = 13.3c 

E 

= 10.35 

E : 

= 10 . 2 lv 

E 

= 13.87 

station 


^ = 

: 180 ° 

^ = 

225 ° 

Mr = 

= 270° 

Mr = 

= 315° 

Fd'ni 

Fd'Hl X 5 " 

Fd'ni 

Fd'Ri X 5 " 

Fd’ni 

Fd'rii X 5 " 

Fd'Hl 

Fd’Hi X 5 " 

0 

-0.210 

0.172 

0.860 

0.111 

0.555 

0.029 

0.145 

-0.026 

-0.130 

1 

-.560 

1.54 

7.70 

1.55 

7.75 

1.07 

5.35 

. 39 *^ 

1.970 

2 

-.680 

2.82 

Ill-.IO 

2.91 

14.55 

1.95 

9.75 

.**97 

2.49 

3 











-.462 

2.32 

11.60 

“2.32 

11.60 

1.33 

6.65 

1 

b 

o^ 

-.325 

Iv 










c 

.010 

-.054 

-.270 

-.051 

-.255 

-.020 

-.100 

.018 

.090 

0 

.610 

-3.13 

-15.65 

-2.59 

-12.95 

-.181 

-.905 

2.72 

13.60 

6 



Z = 18.26 

L 

= 21.25 

E 

1 

= 20.93 

E 

I 

= 17.68 

I 
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TABLE V 


CALCULATION OF 



Station 

’ll 


m* 

( slugs/ft) 

X 5/12 

0 

-0.210 

0.0441 

0.00602 

0.111 X 10“3 

1 

-.560 

.3136 

.00602 

.787 

2 

-.680 

.4624 

.00602 

1.16 

3 

-.462 

.2134 

.00602 

.536 

1 ; 

.010 

.0001 

.00602 

.00025 

5 

.610 

.3721 

.00602 

.934 

6 



pR 





/ m'Ti .2 dx = 3.528 X 10 -^ 




7° J 

1 
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TABLE VI 


CALCULATIONS OF AT EIGHT AZIMUTH POSITIONS 


(deg) 

jTpd’Tii dx 

r ^ 

/ Fd’m dx 
J 0 

.E 2 

/ dx 

J 0 

0 

13.30 

3169 


10.35 

2935 

90 

10. 2U 

2901 

135 

13.87 

3932 

180 

18.26 

5177 

225 

21.25 

602L 

270 

20.93 

5929 

315 

17.68 

5 OIL 


TABIE VII 


CALCULATION OF Si^n fa ^ % 

U> 


= 55 cps; = 62.5 cps; (u = 800 rpm; = 7029. l( radians/sec )^J 


( 1 ) 


( 2 ) 


( 3 ) 


W 


( 5 ) 


( 6 ) 


( 7 ) 


(8) 




®l,nb 


n 2 a^ 


(n,«)" - 


( 5 ) 


gl,na ^ (^) 


81, nb 


0 

1 

2 

3 


4500 

-720 

0 

0 


-1560 

0 

0 


0 

T029.1 


154,056 

147,027 


0.775 

.812 


3490 

-586.0 

0 

0 


-1265 

0 

0 


(1) 

(2) 

( 3 ) 

( 4 ) 

( 5 ) 


d| 

Steady state 

First cosine 

First sine 

i 

%,0 = (’^^n=0 ^ (2) 

Ml,a = (T)n=i X (2) 

Ml,b = (8)n=l X (2) 

0 

0 

0 

0 

0 

1 

.6666 X 10"3 

2.33 

-.391 

-.843 

2 

1.27516 

4.45 

-.75 

-1.61 

3 

1.70116 

5.94 

-.998 

-2.15 

4 

1.89016 

6. 60 

-1.107 

- 2.39 

5 

1.82659 

6.37 

-1.07 

-2.31 

6 

1.53819 

5.37 

-.90 

-1.94 

7 

1.0946 

3.82 

-.64 

-1.38 

8 

.60099 

2.10 

-.35 

-.76 

9 

.1878 

.66 

-.11 

-.24 

10 

0 

0 

0 

0 


'Si 
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TABLE VIII 


ON 


CALCULATION OF 


g2,n X ^ ^ 


{^2y0 cps; V 2 ^(x) = 178.5 cps; o) = 8 OO rpm; od? = 7029.^ 


( 1 ) 

(2) 

( 3 ) 

(1^) 

( 5 ) 

( 6 ) 

(7) 

( 8 ) 

n 

^2,na 

82, nb 


- m 

(v2,o)2 

( 5 ) 

82, na X (6) 

82, nb X (6) 

0 

1 

2 

3 

90 

-947 

0 

0 

- 7^7 

0 

0 

0 

7029.1 

1,257,396.7 

1,250,367.6 

0.8967 

.90177 

80.7 

-853.9 

0 

0 

-673.6 

0 

0 







(1) 

(2) 

( 3 ) 

( 4 ) 

( 5 ) 

i 

JJ^ m*T)2 dx d 5 

Steady state 

“2,0 = (T)n=0 X (2) 

First cosine 
“2,a = ('?)n=l X (2) 

First sine 
“2,b = (S)n=l X (2) 

.. 0 
1 

2-1 

3 

-4 

5 ^ 

6 

7 

8 
9 

10 

0 

-.3818 X 10"3 

-.1*832 
-.5107 
-.1957 
.2166 
i . 5381 * 

.6315 

.1*71*3 

.1878 

0 

0 

-.0308 

-.0389 

-.01*12 

-.0158 

.0175 

. 01*35 

.0509 

,0383 

.0152 

0 

0 

.326 
.413 
. 436 
.167 
-.185 
-.459 

-.539 

-.405 

-.160 

0 

0 

.257 

.325 

.344 

.132 

-.146 

-.363 

-.425 

-.319 

-.126 

0 
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TABLE IX 

SUMMARY OF BENDING-MOMENT DISTRIBUTIONS 
(INCLUDES FIRST AND SECOND NATURAL MODES) 


i 

Mq 

( steady 
state) 

Ma 

( cosine) 
( in. -lb) 

Mb 

( sine) 

( in. -lb) 

Mfflax 

i = 225° 

( in. -lb) 

0 

0 

0 

0 

0 

1 

2.30 

-.06 

-.58 

2.75 

2 

4 Ai 

-.34 

-1.28 

5.55 

3 

5.90 

-.56 

-1.81 

7.58 

k 

6.58 

-.94 

-2.26 

8.84 

5 

6.39 

-1.26 

-2.45 

9.01 

6 

5 . 4 i 

-1.36 

-2.30 

8.00 

7 

3.87 

-1.18 

-1.81 

5.98 

8 

2.14 

-.76 

-1.08 

3.44 

9 

.68 

-.27 

-.37 

1.13 

10 

0 

0 

0 

0 
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TABIE X 

CALCUUTION OF SaEADY-STATE COMPOHENT OF BENDING MOMENT 
FOR A ROTOR BIAIE FIXED AT ROOT 


[u = 0.47; fi = 600 rpm; mi = O.OO602 slug/ft; p = 0 °; Rq = 8°; a = -5°; 
X * 5 in.; El = l6,000 lb-in. Zjjj terms neglected, ( )* indicates 
value from the preceding station (i^ 


i 

(1) 

(2) 

(3) 


(5) 

(6) 

(7) 

^i 

(1)* + (6)» + (7)* 

x2 

Eli 

x2aj 

El 



F 

^■"5 j=i 


0 

0 

1.563 X 10-3 

0 

0 

73.80 

0 

-15.22 

1 

-15.22000 

1.563 

-.0237886 

-.0237886 

69.67 

-1.65736988 

-15.10 

2 

-31.97736 

1.563 

-.0499806 

-.0737692 

61. 4i 

-4.53016836 

-lk .^6 

3 

-51.06753 

1.563 

-.0798185 

-.1535877 

49.02 

-7.5288735 

-13.00 

k 

-71.59641 

1.563 

-.1119052 

-.2654929 

32.50 

-8.628522 

-9.60 

5 

-89.92493 

1.563 

-.1405526 

- . 4060456 

11.85 

-4.8 ii 64094 

-3.70 

6 

-98.43657 

1.563 














(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

1 

t>i 

X^b^ 


F 1(10) 

b^Z 

Mi 


(8)* + (11)* 

Eli 


2 


4- bj^Z 

0 

1 

0.0015630 

0.0015630 

0,1153494 

37.61965 

37.61965 

1 

1.1153494 

.0017432 

.0033062 

.2303493 

41,95906 

26.73906 

2 

1. 3456987 

.0021033 

.0054096 

.3322045 

50.62472 

18.64735 

3 

1.6779033 

.0026225 

.0080321 

.3937376 

63.12214 

12,05467 

4 

2.0716409 

.0032379 

.0112701 

. 3662801 

77.93441 

6, 33800 

5 

2.4379210 

.0038104 

.0150806 

, 1787054 

91.71374 

1,78881 

6 

2.6166264 




98.43657 

0 


2 = - ^ = - -98- ^36^73 ^ 37,619651 
t>6 2.616626 
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TABIE XI 

COMPUTATION CHECK FOR STEADY-STATE COMPONENT OF BENDING MOMENT 
FOR A ROTOR BLADE FIXED AT ROOT 


= O.^Tj n = 600 rpm; nij^ = 0.00602 slug/ft; p = 0°; 6q = 8°; 
a = X = 5 in.j El = 16^000 lb-in. Zjjj terms neglectec^ 



(13) 

(14) 

(15) 

(16) 

(17) 

i 

Mi 

®i + 


X2 

Eli 

x2 

Eli ^ • 


0 

37.61965 

-10.88059 

1.563 X 10-3 

0.0587995 

0.0587995 

1 

26.73906 

-8.09171 

1.563 

.0417932 

. 1005927 

2 

18.64735 

-6.59275 

1.563 

.0291458 

.129739 

3 

12.05467 

-5.71667 

1.563 

.0188413 

. 148580 

k 

6.33800 

-4.54919 

1.563 

. 0099063 

.158486 

5 

1.78881 

-1.78881 

1.563 

.0027959 

. 161282 

6 

0 


1.563 

0 



(18) 

(19) 

(20) 

(21) 

(22) 

i 

L.i 

it| VlfA + i 


m 

(19) + (20) 

yi = ^ (I7)j 

0 

73.80 

4.339991 

-15.22 

-10.88097 

0 

1 

69.67 

7.008293 

-15.10 

-8.09171 

.0587995 

2 

61. 4 i 

7.967272 

-14.56 

-6.59273 

.1593922 

3 

49.02 

7.283391 

-13.00 

-5.71661 

. 2891312 

k 

32.50 

5.150795 

-9.60 

-4.4492 

. 4377112 

5 

6 

11.85 

1.911192 

-3.70 

-1.78881 

.59619(2 

.7574792 


TABLE XII 


VJI 

O 


GOODYEAR METHOD (MODIFIED) FOR CALCULATING FIRST-HARMONIC SINE COMPONENT 

OF BENDING MOMENT FOR A ROTOR BLAIE FIXED AT ROOT 

Qi = OAT; = 600 rpm; mi = 0,00602 slug/ft; 3 = 0°; 0 q = 8°; a = -5°; X = 5 in.; 
El = 16,000 lb-in. ( )* indicates value from preceding station (i^ 

(a) Calculation of a. 


i 

(1) 

(2) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

(8) 

(9) 

(10) 

(11) 

(!)♦ + (6)« + (7)* + (11)* 

x! 

El 

X^ai 

El 

^ ■ J^l 



XS T 

Ui 

(4)* + (8)* 


x^n^iUi 

t 

(10) + (11)* 

0 

0 

1.563 X 10"3 

0 

0 

73.80 

0 

-25.01 

0 

4.13 

0 

0 

1 

-25.01 

1.563 

-.0390906 

-.0390906 

69.67 

-2.7234421 

-25.00 

0 

4.13 

0 

0 

2 

-52.733‘^‘^2l 

1.563 

-.0824224 

-.1215129 

61. 4l 

-7.4621115 

-22.97 

-.0390906 

4.13 

-.1614442 

-.1614442 

3 

-83.3269978 

1.563 

-.1302401 

-.2517530 

49,02 

-12.3409355 

-18.85 

-.1606036 

4.13 

-.6632929 

-.8247371 

k 

-115. 31*26704 

1.563 

-.1802806 

-.4320336 

32.50 

-i4.o4i094o 

-12.75 

-.4123567 

4.13 

-1.7030332 

-2.5277703 

5 

-144.6615347 

1.563 

-.2261060 

-.6581396 

11.85 

-7.7989542 

-4,55 

• -.8443904 

4.13 

-3.4873324 

-6.0151027 

6 

-163.0255917 

1.563 

-.25W090 

-.9129^6 

0 

0 

— 

-1.502530 

4.13 

-6.2054489 

-12.2205516 


(b) Calculation of b. 


1 

i 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(12)*+ (15)* + 
(18)* 

X^bi 

El 

— x2b. 

"‘■Hi/ 

2 

Vi 

(14) + (16) 


(17) + (18)* 

0 

1.0 

0.001563 

0.001563 

0.1153494 

0 

0 

0 

1 

1.1153494 

.001743 

.003306 

,2303499 

.001563 

.0064552 

.0064552 

2 

1.3521545 

.002113 

.005419 

.3328250 

.00^4869 

.0201102 

.0265654 

3 

1.7115449 

.002675 

.008094 

. 3968100 

.010289 

.0424937 

.0690591 

4 

2.1774 i 4 o 

.003403 

.011498 

.3736902 

.018383 

.0759254 

.1449845 

5 

2.6960887 

.004213 

.015712 

^ 1861889 

.029882 

.1234128 

.2683973 

6 

3.1506750 

.004924 

.020636 

0 

.045594 

.1883040 

.4567013 
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TABIE XII - Concluded 


> 


GOODYEAR METHOD (MODIFIED) FOR CALCULATING FIRST-HARMONIC SINE COMPONENT 
OF BENDING MOMENT FOR A ROTOR BIADE FIXED AT ROOT 

(c) Calculation of c. 


1 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

( 24 ) 

(25) 

1.0 + (19)* + (22)* + (25)* 

x 2 ci 

El 

X^c . 

= 2 _^ 
j=l 


Wi 

( 21 )* + (23)* 

x^n^iWi 

± (24)J 

i=l 

( 24 ) + (25)* 

0 

0 

0 

0 

0 

0 

0 

0 

1 

1.0 

. 001563 

.001563 

. 10889421 

0 

0 

1 0 

2 

2. 10889^1 

.003296 

.004859 

.29840347 

.001563 

.0064552 

1 .0064552 

3 

3.1H375288 

.005335 

.010194 

.49975380 

.006422 

.0265236 

i .0329788 

hr 

1 ^. 9461)8557 

.007731 

. 017926 

.58260313 

.016617 

.0686286 

. 1016075 

5 

6.63069620 

.010363 

. 028290 

. 33523686 

.034543 

.1426640 

1 .2442715 

6 

8.21020460 

.012832 

.041122 

0 

.062833 

.2595018 

1 .5037734 


(d) Calculation of M. 


i 

( 26 ) 

(27) 

(28) j 

biMo 

=10 

Mi 

(1) + (26) + ( 27 ) 

0 

147.5 

0 

147.5 

. 1 

164 . 5 i 4 

- 36.8 

102.704 

2 

199.443 

-77.6072 

69.102 1 

3 

252.453 

-125.6261 

43.499 

hr 

321.168 

-182.0306 

23.794 

5 

397.672 

-244.0096 

9.001 

6 

464.725 

-302.1355 

0 



\J1 


TN 2626 


52 


NACA TN 2626 



STAGE 1 - ROLLED AI-(;MIN(,M TIJOE AM) E.OOT INSERT FITTING 



TOP ‘.IDFACK 


P>OTTOM SERFAGE 



STAGE 2 - Sm\IN GACF.S ATI ACHED TO SPAR AM) WIRED 




' . • - 1 

- '• . t 

j 


STACK i - COMPLKTKD BI ABE 


Figure 1.- Blade construction details. 







OF RCr^^T\OH SCI>.UE. 
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R 

1 

2&% 
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<o 
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0.52S 

2. 
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7 

13% 
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3 

Z5^ 
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4 
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Figure 2.- Strain-gage locations on rotor blade. 
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Tlieory (Goodyear Method) 
Experiment 


No Tip Loss (Uniform Infloiw) 


X Variation 
(No Tip Loss) 


^CA 






(a) Hovering, hinged- at- root condition; rotor speed, l400 rpm. 


Figure 6.- Comparison of theoretical and experimental hending-moment 
distributions. 6 = 8^; \ = -0.05^; El = l6,000 pound- inches^. 


3% Tip Loss 
(Uniform Inflow) 
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iTlTTri ITTiTri 


\ Variation 
(No Tip Loss) 


No Tip Loss 
(Unifonn Inflow) 


3% Tip Loss 
(Ihiform Inflow) 




J^ACA 


Wimwm 


iHbl 


(b) Hovering, fixed-at-root condition; 3 q = 0°^ rotor speed, 800 rpm. 


\Sl 

CO 


Figure 6.- Concluded. 
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(a) li = 0.10. 


Figure 7.- Bending-moment traces from oscillograph records. Hinged-at- 
root condition; Oq = 8°; a = -5°; rotor speed, 8 OO rpm. 


6o 
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(t) n = 0.22. 

Figure 7 .- Continued 


Moment moment 
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(c) \i = 0.30. 


Figure 7 *- Continued. 


Moment (.\u.\b.'\ Moment 


62 


NACA TN 2626 



(d) = 0.40. 

41 

Figure 7.- Continued 
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(e) la = 0.50. 


Figure 7 .- Concluded 



Figure 8.- Maximum-bending-moment distribution. Hinged- at -root 
condition; 9 = 8°; a = -5°; t = 235°; rotor speed, 800 rpm; 

p 

El = 16^000 pound-inches . 
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Bending Moment (in. lb.) 



Steady State Component 
1st Harmonic Sine 
1st Harmonic Cosine 
5 th Harmonic Sine 


NOTE: 5th Harmonic Cosine 

Negligible . 




-rttitrii 




Figure 9.- Harmonic- bending- moment distribution. Hinged-at-root condi- 
tion; n = 0.22; rotor speed, 800 rpm; 0 = 8°; a = -5°; X = -0.0358; 
El = 16,000 pound- inches^. 


ON 
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ON 

ON 



Figure 10.- Comparison of theoretical and experimental maximum-bending- 
moment distributions. Hinged-at-root condition; p = 0.22; rotor 
speed, 800 rpm; 9 = 8^; a = -5^; X = -0.0358; El = l6,000 pound- 

2 

inches . 
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Theory 

Steady State Coraponent 
1st Harmonic Sine 
1st Harmonic Cosine 


State (De Guillenchmidt)|^ 


Steady State (Goodyear) 


Steady State (Experiment) 


§5?St'lst Cosine (De GSuiliendimidt ) 
1st Sine (De Cuillenchnidt 


1st Sine (Goodyear) 


imaS!H!!H!!l 


Figure 11.- Comparison of theoretical and experimental harmonic-bending- 
moment distributions. Hinged-at-root condition; p = 0.22; rotor 
speed, 800 rpm; 0 = 8°; a = -5°; X = -O.O358; El = l6,000 pound- 

p 

inches^. 


On 
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(a) |i = 0.15; rotor speed, 6 OO rpm. 

Figure 12.- Bending-moment traces from oscillograph records. Fixed-at- 

root condition; 6 q = 8°; a = -5°- 
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0.47; rotor speed, 6 OO rpm. 
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Figure 12 .- Continued. 
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Figure 12.- Continued. 
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Figure I 3 .- Sample traces from 
condition; 6 = 8°; a = 


oscillograph records. Fixed-at-root 
- 5 °; station 9, r/R = O.296. 
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Figure 14.- 
tion; 


Maximum-bending- moment distribution. Fixed-at-root condi- 
3 = 0^; Q = 8°; a = -5^; El = l6,000 pound- inches^. 
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Figure 15- - Harmonic-lDending-inoinent distrilDution. Fixed-at-root condi- 
tion; \i = 0.47; rotor speedy 600 i^m; p = 0^; 0 = 8^; a = 

X. = -0.0457; El = 16^000 pound- inches^. 
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Figure l6.- Comparison of theoretical and experimental harmonic-bending- 
moment distributions. Fixed-at-root condition; p = 0.47; rotor speed, 
600 rpm; p = 0°; 9 = 8°; a = -5°; X = -0.0457; El = l6,000 pound- 
inches^. 
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Figure 17 - - Variation of first natural bending frequency of rotor blades 
of different stiffnesses with rotational speed. 
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(b) Fixed-at-root condition. 
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(a) Hinged-at-root condition; rotor speed, 1000 rpm. 
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(b) Fixed-at-root condition; rotor speed, 400 rpm. 


Figure l8.- San5)le oscillograph records showing resonance phenomenon. 
9 = 12°; El = 16,000 pound- inches^. 



Rotational Speed (r.p.m.) 


(a) El = 3000 pound- inches^. 


Figure 19.- Rotor-blade resonance study. Hovering, hinged-at-root 

condition; 9 = 12 °. 
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Figure I9.- Continued. 
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Figure 19- - Concluded. 






Figure 20.- Rotor-blade resonance study. Hovering, fixed-at-root 

condition; 6 = 12°. 
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(t) El = 16^000 pound- Inches^. 
Figure 20 .- Continued. 





Figure 20.- Concluded. 
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Figure 21 .- Comparison between rotating and nonrotating natural bending 
modes of the full-scale R-6 rotor blade. 
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Figure 22.- Plot of total external loading on rotor blade at ei^t 
azimuth positions for the De Guillenchmidt method. 
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Figure 23.- Natural modes of rotor blade. Nonrotating; El = 16^000 poimd- 

inches^; m* = O.OO 603 slug per foot. 
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Figure 24.- Curves of g-, and gp for De Guillenchmidt method. 
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NOTE: Numbers Oh Curves 
Indicate the Stage 
of Iteration. 
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Figure 25 .- Plot showing convergence of modified Goodyear method for 
calculating bending moment on a rigid rotor blade. 
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